Statistics of Impedance and Scattering Matrices 
in Chaotic Microwave Cavities: Single Channel Case 
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We study the statistical properties of the impedance (Z) and scattering (S) matrices of open 
electromagnetic cavities with several transmission lines or waveguides connected to the cavity. In 
this paper, we mainly discuss the single port case. The generalization to multiple ports is treated in 
a companion paper. The model we consider is based on assumed properties of chaotic eigenfunctions 
for the closed system. Analysis of the model successfully reproduces features of the random matrix 
model believed to be universal, while at the same time incorporating features which are specific 
to individual systems as treated by the Poisson kernel of Mello et al. Statistical properties of the 
cavity impedance Z are obtained in terms of the radiation impedance (i.e., the impedance seen at a 
port with the cavity walls moved to infinity). Efi^ects of wall absorption are discussed. Theoretical 
predictions are tested by direct comparison with numerical solutions for a specific system. (Here 
the word universal is used to denote high frequency statistical properties that are shared by the 
members of the general class of systems whose corresponding ray trajectories are chaotic. These 
universal properties are, by definition, independent of system-specific details.) 
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I. INTRODUCTION 

The problem of the coupling of electromagnetic radia- 
tion in and out of structures is a general one which finds 
applications in a variety of scientific and engineering con- 
texts. Examples include the susceptibility of circuits to 
electromagnetic interference, the confinement of radia- 
tion to enclosures, as well as the coupling of radiation to 
structures used to accelerate charged particles. 

Because of the wave nature of radiation, the coupling 
properties of a structure depend in detail on the size and 
shape of the structure, as well as the frequency of the 
radiation. In considerations of irregularly shaped elec- 
tromagnetic enclosures for which the wavelength is fairly 
small compared with the size of the enclosure, it is typical 
that the electromagnetic field pattern within the enclo- 
sure, as well as the response to external inputs, can be 
very sensitive to small changes in frequency and to small 
changes in the configuration. Thus, knowledge of the re- 
sponse of one configuration of the enclosure may not be 
useful in predicting that of a nearly identical enclosure. 
This motivates a statistical approach to the electromag- 
netic problem. 

While our ability to numerically compute the response 
of particular structures has advanced greatly in recent 
years, the kind of information needed for a statistical de- 
scription may not be obtainable directly from numerical 
computation. 

Thus it would seem to be desirable to have specific an- 
alytical predictions for the statistics of electromagnetic 
quantities in such circumstances. This general prob- 
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lem has received much attention in previous work (e.g., 
Refs. H S H S S B 0)- Some of the main issues ad- 
dressed are: the probability distribution of fields at a 
point, the correlation function of fields at two points near 
each other, the statistics of the excitation of currents in 
cables or in small devices within the enclosure, the cav- 
ity Q, the statistics of coupling to the enclosure, and the 
statistics of scattering properties. A fundamental basis 
for most of these studies is that, due to the complexity of 
the enclosure and the smallness of the wavelength com- 
pared to the enclosure size, the electromagnetic fields ap- 
proximately obey a statistical condition that we shall call 
the random plane wave hypothesis, which assumes that a 
superposition of random plane wave can be used to de- 
scribe the statistics of chaotic wave functions 8]. This 
work has been quite successful in obtaining meaningful 
predictions, and some of these have been tested against 
experiments with favorable results. A good introduction 
and overview is provided in the book by Holland and St. 
John 0. 

In addition to this previous work on statistical elec- 
tromagnetics P, S S IJj S S 1 much related work has 
been done by theoretical physicists. The physicists are 
interested in solutions of quantum mechanical wave equa- 
tions when the quantum mechanical wavelength is short 
compared with the size of the object considered. Even 
though the concern is not electromagnetics, the questions 
addressed and the results are directly applicable to wave 
equations, in general, and to electromagnetics, in par- 
ticular. The start of this line of inquiry was a paper 
by Eugene Wigner 0. Wigner's interest was in the en- 
ergy levels of large nuclei. Since the energy level den- 
sity at high energy is rather dense, and since the solu- 
tion of the wave equations for the levels was inaccessible, 
Wigner proposed to ask statistical questions about the 
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levels. Wigner's results apply directly to the statistics 
of resonant frequencies in highly-overmoded irregularly- 
shaped electromagnetic cavities. Since Wigner's work, 
and especially in recent years, the statistical approach to 
wave equations has been a very active area in theoretical 
physics, where the field has been called 'quantum chaos'. 
We emphasize, however, that the quantum aspect to this 
work is not inherent, and that a better terminology, em- 
phasizing the generality of the issues addressed, might be 
'wave chaos'. For a review see Chapter 11 of Ref. [Ty| or 
the books [Hill- 

Wigner's approach was to introduce what is now called 
Random Matrix Theory (RMT) In RMT the hnear 
wave equation is replaced or modelled by a linear matrix 
equation where the elements of the matrix are random 
variables. This follows from Wigner's hypothesis that the 
eigenvalues for a complicated (in our case chaotic) system 
have the same statistics as those of matrices drawn from 
a suitable ensemble. Based on symmetry arguments, 
Wigner proposed that the matrix statistics are those that 
would result if the matrix were drawn from different types 
of ensembles, where the relevant ensemble type depends 
only on gross symmetries of the modelled system. The 
two ensembles that are relevant to electromagnetic prob- 
lems are the Gaussian Orthogonal Ensemble (GOE) and 
the Gaussian Unitary Ensemble (GUE). In both cases, 
all the matrix elements are zero mean Gaussian random 
variables. In the GOE all the diagonal element distri- 
butions have the same width, while all the off diagonal 
element distributions have widths that are half that of 
the diagonal elements. The matrices are constrained to 
be symmetric, but otherwise the elements are statisti- 
cally independent. The GOE case is intended to model 
wave systems that have time reversal symmetry (TRS). 
That is, the time domain equations are invariant under 
the transformation t —t. This is the case for elec- 
tromagnetic waves if the permitivities and permeabilities 
tensors are real and symmetric. In the GUE the ma- 
trices are constrained to be Hermitian. In this case the 
off-diagonal elements are complex and the distributions 
of their real and imaginary parts are independent and 
Gaussian and the width of these Gaussians is again the 
one half the width of the real diagonal elements. The 
GUE case is intended to model systems for which time 
reversal symmetry is broken (TRSB). This case will ap- 
ply in electromagnetics if a nonreciprocal element such 
as a magnetized ferrite or a cold magnetized plasma, is 
added to the system. 

In this paper we mainly consider an irregularly shaped 
cavity with a single transmission line and/or waveguide 
connected to it, and we attempt to obtain the statistical 
properties of the impedance Z and the scattering matrix 
S (which are both scalars in the cases we consider) char- 
acterizing the response of the cavity to excitations from 
the connected transmission line, where the wavelength is 
small compared to the size of the cavity. We will treat 
specifically the case of cavities that are thin in the ver- 
tical (z-direction) direction. In this case the resonant 



fields of the closed cavity are transverse electromagnetic 
{TMz, E = Ez{x, y)z), and the problem admits a purely 
scalar formulation. While the two dimensional problem 
has practical interest in appropriate situations (e.g., the 
high frequency behavior of the power plane of a printed 
circuit), we emphasize that the results for the statistical 
properties of Z and S matrices are predicted to apply 
equally well to three dimensional electromagnetics and 
polarized waves. We note that previous work on sta- 
tistical electromagnetics d SH, MlE IE is for fully 
three dimensional situations. Our main motivation for 
restricting our considerations here to two dimensions is 
that it makes possible direct numerical tests of our pre- 
dictions (such numerical predictions might be problem- 
atic in three dimensions due to limitations on computer 
capabilities). Another benefit is that analytical work and 
notation are simplified. 

For an electrical circuit or electromagnetic cavity with 
ports, the impedance matrix provides a characterization 
of the structure in terms of the linear relation between 
the voltages and currents at all ports, 

V = zi, (1) 

where V and / are column vectors of the complex phasor 
amplitudes of the sinusoidal port voltages and currents. 
The scattering matrix S is related to the impedance ma- 
trix Z by 

S = Z',/\Z + Zo)-\Z - Zo)Z,'^\ (2) 

where Zq is the characteristic impedance of the transmis- 
sion line. 

As discussed in the next section, the impedance ma- 
trix Z can be expressed in terms of the eigenfunctions 
and eigenvalues of the closed cavity. We will argue that 
the elements of the Z matrix can be represented as com- 
binations of random variables with statistics based on the 
random plane wave hypothesis for the representation of 
chaotic wave functions and results from random matrix 
theory for the distribution of the eigenvalues. 

This approach to determination of the statistical prop- 
erties of the Z and S matrices allows one to include the 
generic properties of these matrices, as would be pre- 
dicted by representing the system as a random matrix 
drawn from an approapriate ensemble. It also, however, 
allows one to treat aspects of the S and Z matrices which 
are specific to the problems under consideration (i.e., so- 
called non-universal properties). 

These nonuniversal properties have previously been 
treated within the context of the so-called Poisson Ker- 
nel based on a "maximum information entropy" principle 
'l4l, and Brouwcr later provided a microscopic justifica- 
tion and showed that the Poisson Kernel can be derived 
from Wigner's RMT description of the Hamilton [l5l| . 
Here the statistics of the S matrix depend in a non-trivial 
way on the average of S taken in a narrow frequency 
range. This characterizes the system specific aspects of 
the coupling. Our approach allows one to predict the 
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average based on another informative quantity, the radi- 
ation impedance, which itself characterizes the coupHng 
of the port to the enclosure. The radiation impedance 
is the impedance that applies at the port when waves 
are launched into the cavity and (by making the distant 
walls perfectly absorbing) not allowed to return. Our 
interpretation of the role of the radiation impedance is 
equivalent to Brouwer's (15>] interpretation of the Pois- 
son kernel in terms of scattering from a cascaded con- 
figuration of a lossless multiport and a perfectly coupled 
cavity described by Random Matrix Theory. We show 
that the separation of universal and system specific prop- 
erties is more natural when considering the impedance 
rather than scattering matrix. Specifically, the univer- 
sal properties of the cavity impedance are observed by 
subtracting from the raw cavity impedance the radiation 
reactance and normalizing the result to the radiation re- 
sistance. The statistics of the resulting variable, which 
we term the normalized cavity impedance, depend only 
on a single parameter characterizing the internal loss of 
the cavity. 

The Poisson kernel description of the scattering am- 
plitude has been applied to data obtained from mi- 
crowave scattering experiments on cavities with absorp- 
tion 16, 17, l&l. Two different methods have recently 
been described |l7l Il8l| for extracting the universal prop- 
erties of the scattering amplitude from the system spe- 
cific ones. The system specific properties are character- 
ized in terms of averages of the scattering amplitude and 
its modulus squared over ranges of frequency. A fitting 
procedure is then used to characterize the universal fluc- 
tuations. The impedance approach has also been applied 
to experimental data 20} ■ In I'ef. the radiation 
impedance is measured directly and the complex nor- 
malized impedance is formed by subtracting the radia- 
tion reactance and normalized to the radiation resistance. 
The probability distribution function of the normalized 
impedance was then compared directly with Monte Carlo 
numerical evaluations of the prediction. 

The main contribution of this paper are to describe 
the statistical properties of the cavity impedance. First 
we show (also shown in Ref. that the relation be- 
tween the cavity impedance and the radiation impedance 
follows from the assumption that the eigenfunctions of 
the cavity satisfy the random plane wave hypothesis of 
Berry We then verify the relation between the cav- 
ity impedance and the radiation impedance for a specific 
realization by numerical simulation using the High Fre- 
quency Structure Simulation software, HFSS. Next, using 
Monte Carlo methods we evaluate theoretical predictions 
for the probability distribution functions for the real and 
imaginary parts of the normalized impedance for differ- 
ent values of internal absorption in the cavity. The dis- 
tribution of the real part of the normalized impedance is 
closely related to the distribution of values that is known 
as the local density of states . Finally, we derive ex- 
pressions for the mean and variance of the normalized 
impedance as functions of the level of internal loss in the 



cavity. 

Our paper is organized as follows. In Sec. II, we 
presents the statistical model. Section III illustrates our 
model by application to the statistics of the impedance 
seen at a single transmission line input to a cavity that is 
irregularly-shaped, highly over-moded, lossless, and non- 
gyrotropic (i.e., no magnetized ferrite). Section IV re- 
lates the impedance matrix characteristics to those of 
the scattering matrix. Section V generalizes our model 
to incorporate the effects of distributed loss (such as wall 
absorption). Throughout, our analytical results will be 
compared with direct numerical solutions of the wave 
problem. Section VI concludes with a discussion and 
summary of results. 



II. MODELLING WITH RANDOM PLANE 
WAVES 

We consider a closed cavity with ports connected to 
it. For specificity, in our numerical work, we consider the 
particular, but representative, example of the vertically 
thin cavity shown in Fig. ^a) coupled to the outside via 
a coaxial transmission cable. Fig.^b) shows an example 
of how this cavity might be connected to a transmission 
line via a hole in the bottom plate. The cavity shape in 
Fig. n is of interest here because the concave curvature 
of the walls insures that typical ray trajectories in the 
cavity are chaotic. (Fig.QIa) is a quarter of the billiard 
shown in Fig. 2(c).) For our purposes, a key consequence 
of the chaotic property of the shape in Fig 1(a) is that, 
if we consider the trajectory of a particle bouncing with 
specular reflection at the walls (equivalently a ray path), 
then a randomly chosen initial condition (i. e., random in 
position X within the cavity and isotropically random in 
the orientation 9 of the initial velocity vector) always gen- 
erates an orbit that is ergodic within the cavity. Here by 
ergodic we mean the following: For any spatial region R 
within the cavity, in the limit of time t — > oo, the fraction 
of time the orbit spends in R is the ratio of the area of R 
to the entire area of the two-dimension cavity, and, fur- 
thermore, the collection of velocity orientations of the 
orbit when it is in R generates a uniform distribution in 
[0, 27r). Thus the orbit uniformly covers the phase space 
(x, 6). This is to be contrasted with the case of a rect- 
angular cavity, 0<x<a, 0<y<6, which represents 
a shape for which orbits (rays) are nonchaotic. In that 
case, if the initial velocity orientation with respect to the 
X-axis is 9q then at any subsequent time only four values 
of 9 are possible: 27r — 0o, — ^'oi 7^ + 9^. In cases such 
as Fig. 1(a) we assume that the previously mentioned 
hypotheses regarding eigenfunctions and eigenvalue dis- 
tributions provide a useful basis for deducing the statis- 
tical properties of the Z and S matrices, and, in what 
follows, we investigate and test the consequences of this 
assumption. 

The vertical height h of the cavity is small, so that, for 
frequencies of interest, the only propagating waves inside 
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(a) Top View 



a 




R2=63.5cm 



42.5 cm 

(b) Side View of Coupling Details 



h=0.2 cm 




Inner Radius 0. 1 cm 



Outer Radius 0.25 cm 



Coaxial Transmission Line 



FIG. 1: (a) Top view of the cavity used in our numerical 
simulation, (b) Side view of the details of a possible coupling. 



the cavity have electric fields that are purely vertical, 

E^E,{x,y)z. (3) 

This electric field gives rise to a charge density on the 
top plate ps — —eoEz, and also generates a voltage 
Vt{x, y) — —hEz{x, y) between the plates. The magnetic 
field is perpendicular to £, 



B — {Bx, By) — fioH, 



(4) 



and is associated with a surface current density Jg = 
H X z flowing on the top plate. 

The cavity excitation problem for a geometry like that 
in Fig. ^b) is system specific. We will be interested 
in separating out statistical properties that are inde- 
pendent of the coupling geometry and have a universal 
(i.e., system- independent) character. For this purpose, 
we claim that it suffices to consider a simple solvable 
excitation problem, and then generalize to more compli- 
cated cases, such as the coupling geometry in Fig. ^b). 
Thus we consider the closed cavity (i.e., with no losses 
or added metal), with localized current sources result- 
ing in a current density Js{x,y,t) — '^^Ii{t)ui{x,y)z 
between the plates. The profile functions Ui{x, y) are as- 
sumed to be localized; i.e., Ui(x,y) is essentially zero for 
{x — Xif + {y — yiY > If, where k is much smaller than 
the lateral cavity dimension. Ui(x,y) characterizes the 



distribution of vertical current at the location of the i- 
th model input (analogous to the i-th transmission line 
connected to the cavity, although, for this model there 
are no holes in the upper or lower plates). The profile is 
normalized such that 



dxdyui{x, y) = 1. 



(5) 



For the sake of simplicity, we only consider the single 
port case in this paper (i.e., there is only one localized 
source and we may thus drop the subscript i on Ui(x, y)). 
The injection of current serves as a source in the continu- 
ity equation for surface charge, dps/dt+V ■ Jg — Iu{x, y), 
where V — {d/dx,d/dy). Expressed in terms of fields, 
the continuity equation becomes: 



d_ 

di 



i^eoE,) + V ■ (H X z) ^ Iu{x,y). 



(6) 



Differentiating Eq. Q with respect to t and using Fara- 
day's law, we obtain, 

^{-eoEz) + V • -VEz = u{x,y)^. (7) 

Expressing the electric field in terms of the voltage Vr = 
—Ezh, we arrive at the driven wave equation. 



1^ 



— V Vt ~ hpou 



5/ 



(8) 



where c is speed of light, = l/(/zoeo). 

Assuming sinusoidal time dependence e^'^* for all field 
quantities, we obtain the following equation relating Vt 
and /, the phasor amplitudes of the voltage between the 
plates and the port current. 



(V^ 4- /c^)Vt = —jujhpoui = —jkhrjoul, 



(9) 



where t/q = \/po/eQ is the characteristic impedance of 
free space and k = uj/c. Thus Eq. © represents a wave 
equation for the voltage between the plates excited by 
the input current. 

To complete our description and arrive at an expression 
of the form of Eq. Q, we need to determine the port 
voltage V. We take its definition to be a weighted average 
of the spatially dependent voltage VT{x,y,t), 



V= dxdyu{x,y)VT{x,y,t). 



(10) 



This definition is chosen because it then follows from 
Eq. that the product IV gives the rate of change 
of field energy in the cavity, and thus Eq. ((Tn|l provides a 
reasonable definition of port voltage. Solution of Eq. Q 
and application of IjlOl) to the complex phasor amplitude 
Vt provide a linear relation between V and /, which de- 
fines the impedance Z . 

To solve Eq. we expand Vr in the basis of the 
eigenfunctions of the closed cavity, i.e., Vt = J^n'^ri'Pn, 
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where (V^ + fc^)(/!)„ = 0, J (f>i(j)jdxdy = Sij and (j)n{x,y) = 
at the cavity boundary. Thus, multiplying Eq. @ by 
(f>n and integrating over (x, y) yields 



(11) 



where /c„ = i-^n/c, ujn is the eigenfrequency associated 
with 0„, and ( u(f>n) = J (finudxdy. Solving for the coef- 
ficients c„ and computing the voltage V yields 



Y _ _ . khi]o{u(j)n)'^ 



ZI. 



(12) 



will be a Gaussian random variable with zero mean. The 
variance of the overlap integral can be obtained using 
Eq. 

E{{u^,,r}^^l^^ ^H^r, (14) 

where E{.} denotes the expected value, u(fc„) is the 
Fourier transform of the profile function u{x,y), 



u{kn) ^ dxdyu{x,y)exp{-jk„ ■ x), (15) 



This equation describes the linear relation between the 
port voltage and the current flowing into the port. Since 
we have assumed no energy dissipation so far (e.g., due 
to wall absorption or radiation), the impedance of the 
cavity is purely imaginary, as is indicated by Eq. (|12|l . 

The expression for Z in Eq. H12|l is equivalent to a 
formulation introduced by Wigner and Eisenbud 22] in 
nuclear-reaction theory in 1947, which was generalized 
and reviewed by Lane and Thomas , and Mahaux and 
Weidenmiiller |24| . Recently, a supersymmetry approach 
to scattering based on this formulation was introduced 
by Verbaarschot et.al. and further developped by 

Lewenkopf and Fyodorov 27] (which they called the 
"if-matrix" formalism), and it has also been adapted to 
quantum dots by Jalabert, Stone and Alhassid |2a ]. 

Explicit evaluation of Eq. H12|) in principle requires de- 
termination of the eigenvalues and corresponding eigen- 
functions of the closed cavity. We do not propose to 
do this. Rather, we adopt a statistical approach to re- 
place (u0n) and fc^ with random variables with appro- 
priate distribution, such that we can construct models 
for the statistical behavior of the impedance. For high 
frequencies such that k = lu/c ^ where L is a typ- 
ical dimension of the cavity, the sum in Eq. H12|l will 
be dominated by high order (short wavelength) modes 
with fc„L ^ 1, and the properties of the short wave- 
length eigenfunctions can be understood in terms of ray 
trajectories. For geometries like that in Fig. ^a), ray 
trajectories are chaotic. 

The assumed form of the eigenfunction from the ran- 
dom plane wave hypothesis is 



lim 

AT— >oo 



N 



AN 



i?e{^ Ui exp(jfc„ei • x + j9i)}, (13) 



i=l 



where Ci are randomly oriented unit vectors (in the x-y 
plane), 9i is random in [0, 27r], and ai are random. This 
statistical model for 0„ is motivated by the previously 
discussed ergodicity of ray paths in chaotic cavities (e.g.. 
Fig. 1(a)); i.e., the random orientation of corresponds 
to the uniform distribution of ray orientations 9. Using 
(|13|l we can calculate the overlap integral {u(j)n} appear- 
ing in the numerator of H12|) . Being the sum of contri- 
butions from a large number of random plane waves, the 
central limit theorem implies that the overlap integral 



and kn = (fc„ cos 6*, fc„ sin 6*). The integral in H14() over 9 
represents averaging over the directions ej of the plane 
waves. The variance of {u(f>n) depends on the eigenvalue 
k^. If we consider a localized source u{x,y) such that 
the size of the source is less than the typical wavelength 
2n/kn, then the variance will be (recall the normal- 
ization of u given by Eq. lO). 

Modelling of Eq. H12|l also requires specifying the dis- 
tribution of eigenvalues A:„ appearing in the denominator. 
According to the Weyl formula 10], for a two dimensional 
cavity of area the average separation between adjacent 
eigenvalues, fc^ — fc^_]^, is AttA'^. Thus, one requirement 
on the sequence of eigenvalues is that they have a mean 
spacing AttA^^. The distribution of spacings of adjacent 
eigenvalues is predicted to have the characteristic Wigner 
form for cavities with chaotic trajectories. In particular, 
defining the normalized spacing, s„ = ^(^n ~ ^^_i)/47r, 
it is found that there are two basic cases which (for rea- 
sons explained subsequently) are called "time reversal 
symmetric" (TRS) and "time-reversal symmetry broken" 
(TRSB). The probability density function for s„ is pre- 
dicted to be closely approximated by 



P{sn) = -s„exp(-7r4/4) 



(16) 



for chaotic systems with time-reversal symmetry(TRS) 
and 



32 

-P(s„) = — 4exp(-4s2/7r) 

TT 



(17) 



for time-reversal symmetry broken(TRSB) system. 
Thus, a second requirement on the sequence of eigen- 
values is that they have the correct spacing distribution. 
The TRS case applies to systems where the permittiv- 
ity and permeability tensors are real and diagonal. The 
TRSB case applies to systems where the permittivity or 
permeability tensors are complex but hermitian, as they 
are for a magnetized ferrite. 

One approach of ours will be to generate values for 
the impedance assuming that sequences of eigenvalues 
can be generated from a set of separations s„ which are 
independent and distributed according to Eq. I|16() . The 
usefulness of the assumption of the independence of sepa- 
rations will have to be tested, as it is known that there are 
long range correlations in the spectrum, even if nearby 
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eigenvalues appear to have independent spacings. A more 
complete approach is to use a sequence of eigenvalues 
taken from the spectra of random matrices. When this 
is done the impedance defined in Eq. 112|) (with indepen- 
dent Gaussian distributions for the overlap integrals) is 
completely equivalent to that obtained in Random Ma- 
trix Theory. We will find that in some cases it is sufficient 
to consider the simpler spectra, generated from indepen- 
dent spacing distributions, but in other cases, for exam- 
ple, when losses are considered, or when correlations of 
impedance values at different frequencies are considered, 
the correlations in eigenvalues exhibited by random ma- 
trix theory are important. This will be discussed more 
thoroughly later in the paper. 

A key assumption in our model is the statistical inde- 
pendence of the overlap integrals, {u(j)n), and the eigen- 
values fc„. This we argue on the basis that each eigen- 
function satisfies the plane wave hypothesis and succes- 
sive eigenf unctions appear to be independent. A second 
justification comes from random matrix theory where it is 
known that the probability distribution for the eigenval- 
ues of a random matrix is independent of that of the ele- 
ments of the eigenfunctions (|13|, Chap. 3). Indeed, the 
result from the random plane wave hypothesis (Eq. (fH^ . 
below) turns out to be equivalent to past work on scat- 
tering matrices that was based on coupling to systems 
described by random matrix theory jl4| . 

Combining our expressions for (u0„) and using the re- 
sult that for a two dimensional cavity the mean spacing 
between adjacent eigenvalues is A = AttA'^, the expres- 
sion for the cavity impedance given in Eq. (|12() can be 
rewritten, 

n—l 

where w„ is taken to be a Gaussian random variable with 
zero mean and unit variance, the fc„ are distributed in- 
dependent of the Wn, and Rr is given by 

R-i^)-'-^jfj<^)\'- (19) 

Our rationale for expressing the impedance in the form 
of Eq. H18() and introducing RR{kn) is motivated by 
the following observation. Suppose we allow the lateral 
boundaries of the cavity to be moved infinitely far from 
the port. That is, we consider the port as a 2D free-space 
radiator. In this case, we solve Eq. (|5J| with a boundary 
condition corresponding to outgoing waves, which can be 
readily done by the introduction of Fourier transforms. 
This allows us to compute the phasor port voltage V by 
Eq. H10|l . Introducing a complex radiation impedance 
Znik) = V/I (for the problem with the lateral bound- 
aries removed), we have 



where RR{kn) is given by Eq. H19() and fc„ is now a con- 
tinuous variable. The impedance Zii{k) is complex with 
a real part obtained by deforming the fc„ integration con- 
tour to pass above the pole at fc„ = k. This follows as 
a consequence of applying the outgoing wave boundary 
condition, or equivalently, letting k have a small nega- 
tive imaginary part. Thus, we can identify the quantity 
Rii{k) in Eq. H19|l as the radiation resistance of the port 
resulting from one half the residue of the integral in H20() 
at the pole, fc^ = 

Re[ZR{k)] = Rnik), (21) 

and 

XR{k) = Im[ZR{k)] 

is the radiation reactance given by the principal part (de- 
noted by P) of the integral (^0)) . 

XR[k)^P{--j^ ^^^R^ik^)}. (22) 

Based on the above, the connection between the cav- 
ity impedance, represented by the sum in Eq. H18|l . 
and the radiation impedance, represented in Eq. H21(l 
and Eq. I|22|l . is as follows. The cavity impedance, 
Eq. (|18|1 . consists of a discrete sum over eigenvalues fc„ 
with weighting coefficients Wn which are Gaussian ran- 
dom variables. There is an additional weighting factor 
Ru{kn) in the sum, which is the radiation resistance. 
The radiation reactance, Eq. (|22|l . has a form analogous 
to the cavity impedance. It is the principle part of a con- 
tinous integral over /c„ with random coupling weights set 
to unity. While, Eqs. ijlHl), (jHl), have been obtained 
for the simple model input J = /u(a;, y) in < z < /i with 
perfectly conducting plane surfaces at z = 0, h, we claim 
that these results apply in general. That is, for a case 
like that in Fig. 1(b), Zji{k) (which for the simple model 
is given by Eq. (HOJl) can be replaced by the radiation 
impedance for the problem with the same port geometry. 
We also note that while H2U|) was obtained with refer- 
ence to a two dimensional problem, the derivation and 
result are the same in three dimensions. It is important 
to note that, while Rn^k) is nonuniversal (i.e., depends 
on the specific coupling geometry, such as that in Fig. 
2(b)), it is sometimes possible to independently calculate 
it, and it is also a quantity that can be directly measured 
(e.g., an experimental radiation condition can be simu- 
lated by placing absorber adjacent to the lateral walls). 
In the next section, we will use the radiation impedance 
to normalize the cavity impedance yielding a universal 
distribution for the impedance of a chaotic cavity. 

III. IMPEDANCE STATISTICS FOR A 
LOSSLESS, TIME REVERSAL SYMMETRIC 
CAVITY 

In the lossless case, the impedance of the cavity Z in 
Eq. (|18|l is a purely imaginary number and S, the reflec- 
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tion coefScient, is a complex number with unit modulus. 
Terms in the summation of Eq. H18|) for which is close 
to k'^ will give rise to large fluctuations in Z as either fc^ 
is varied or as one considers different realizations of the 
random numbers. The terms for which is far from k'^ 
will contribute to a mean value of Z. Accordingly, we 
write 



z ^ z + z, 

where Z, the mean value of Z , is written as 



IE 



AE{ 



1.2 1.2 

1% rv^-, 



(23) 



(24) 



and we have used the fact that the w'^ are independent 
with E{w^} = 1. If we approximate the summation in 
Eq. 124|) by an integral, noting that A is the mean spacing 
between eigenvalues, comparison with H22I) yields 



Z^jXnik), 



(25) 



where Xj^ = Im[Zji] is the radiation reactance defined 
by Eq. H22|l . Thus, the mean part of the fluctuating 
impedance of a closed cavity is equal to the radiation 
reactance that would be obtained under the same cou- 
pling conditions for an antenna radiating freely; i.e., in 
the absence of multiple reflections of waves from the lat- 
eral boundaries of the cavity. The equivalent conclusion 
for the radiation scattering coefRcient is evident from the 
treatment of Brouwer . 

We now argue that, if fc^ is large enough that many 
terms in the sum defining Z satisfy fc^ < fc^, then the 
fluctuating part of the impedance Z has a Lorentzian 
distribution with a characteristic width Rji{k). That is, 
we can write 



Z^jiXR + RiiO. 



(26) 



where ^ is a zero mean unit width Lorentzian distributed 
random variable, P^{£,) = [7r(I +^^)]^^. 

Lorentzian distribution appears in the theory of nu- 
clear scattering .301 and arises as consequences of random 
matrix theory [23, That the characteristic width 

scales as RR{k) follows from the fact that the fluctuating 
part of the impedance is dominated by terms for which 
fc^ ~ k^. The size of the contribution of a term in the 
sum in Eq. IjlHI) decreases as — in the denominator 
increases. The many terms with large values of |fc^ — fc^| 
contribute mainly to the mean part of the reactance with 
the fluctuations in these terms cancelling one another due 
to the large number of such terms. The contributions to 
the mean part from the relatively fewer terms with small 
values of — fc^l tend to cancel due to the sign change 
of the denominator while their contribution to the fluc- 
tuating part of the reactance is significant since there are 
a smaller number of these terms. Consequently, when 
considering impedance fluctuations, it suffices to treat 
RB.{kn) as a constant in the summation in Eq. (|f 8|l and 
factor it out. This results in a sum that is independent 
of coupling geometry and is therefore expected to have a 
universal distribution. 



A. Numerical Results for a Model Normalized 
Impedance 

To test the arguments above, we consider a model nor- 
malized cavity reactance ^ = X/Rji and also introduce a 
normalized wavenumber k'^ — fc^/A — k'^A/Air. In terms 
of this normalized wavenumber, the average of the eigen- 
value spacing [average of {k'^^i — fc^)] is unity. Our model 
normalized reactance is 



N 



1 w„ 



(27) 



where the w„ are independent Gaussian random vari- 
ables, k^ are chosen according to various distributions, 
and we have set i?ij(fc„) to a constant value for n < N 
and R]i{kn) — for n > N. The fluctuating part of 
given by Eq. H27|l mimics the fluctuating part of the 
impedance Z in the case in which i?_R(fc„) has a sharp cut- 
off for eigenmodes with n > N. In terms of ^, Eq. H26fl 
becomes 



1 



1 



(28) 



where f is the mean of ^. 

First we consider the hypothetical case where the col- 
lection of k^ values used in Eq. H27|l result from N in- 
dependent and uniformly distributed random choices in 
the interval ^ fc^ ^ TV. In contrast to Eqs. lfTC|l . 
this corresponds to a Poisson distribution of spacings 
P(s) = exp{—s) for large N. This case is analytically 
solvable j29j and that the mean value ^ is 



1 



N 



1. 

n\—. 1, (29) 



fc2 



and, furthermore, that ^ has a Lorentzian distribution 
given by Eq. if^ . 

Our next step is to numerically determine the proba- 
bility distribution function for ^ given by (|27|l in the case 
where the spacing distribution corresponds to the TRS 
case described by Eq. (3). We generated 10^ realizations 
of the sum in Eq. H27() . For each realization we randomly 
generated = 2000 eigenvalues using the spacing proba- 
bility distribution (3), as well as iV = 2000 random values 
of Wn chosen using a Gaussian distribution for Wn with 
E{wn} = and E{wf^} = 1. We first test the prediction 
of Eq. (|29() by plotting the median value of ^ versus k^ 
in Fig. |2Ia). (We use the median rather than the mean, 
since, for a random variable with a Lorentzian distribu- 
tion, this quantity is more robust when a finite sample 
size is considered.) Also plotted in Fig. EJa) is the for- 
mula (|29|l . We see that the agreement is very good. Next 
we test the prediction for the fluctuations in ^ by plotting 
a histogram of ^ values for the case fc^ = ]V/2 in Fig.[21[b). 
From H29|) for k^ = N/2 the mean is expected to be zero, 
and, as can be seen in the figure, the histogram (open 
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circles) corresponds to a Lorentzian with zero mean and 
unit width (sohd hne) as expected. Histograms plotted 
for other values of fc^ agree with the prediction but are 
not shown. Thus, we find that the statistics of £, are the 
same for P{s) = exp{—s) (Poisson) and for P{s) given 
by Eq. Hence we conclude that the statistics of ^ 

are independent of the distribution of spacings. This is 
further supported by Fig. ^c) where the histogram of f 
for — N/2 is plotted for the case in which the spacing 
distribution is that corresponding to time reversal sym- 
metry broken (TRSB) systems, (the TSRB case will be 
discussed more carefully in a subsequent paper). Again 
the histogram is in excellent agreement with H28|l . This 
implies that, for the lossless case, with a single input 
transmission line to the cavity, the impedance statistics 
are not so sensitive to the spacing distributions, as long 
as they have the same mean value. 

The issue of long range correlations in the distribu- 
tion of eigenvalues seems doesn't affect statistics of the 
impedance in the lossless case. In principle, one can also 
incorporate additional eigenvalue correlation from ran- 
dom matrix theory in the statistics generating the fc^ 
in Eq. H27|l .fand when losses are considered, this is nec- 
essary.) We note that the mean and width of the dis- 
tribution in the random matrix approach are specific to 
the random matrix problem. In contrast, in our formu- 
lation, these quantities are determined by the geometry 
specific port coupling to the cavity through the radiation 
impedance Zfi{k'^). 



B. HFSS simulation result for the normalized 
impedance 

To test our prediction for the distribution function 
of the normalized impedance, we have computed the 
impedance for the cavity in Fig. H^a) for the coupling 
shown in Fig. ^b) using the commercially available pro- 
gram HFSS (High Frequency Structure Simulator [31|). 
To create different realizations of the configuration, we 
placed a small metallic cylinder of radius 0.6 cm and 
height h at 100 different points inside the cavity. In ad- 
dition, for each location of the cylinder, we swept the 
frequency through a 2.0 GHz range (about 100 modes) 
from 6.75GHz to 8.75GHz in 4000 steps of width 5 x lO"'' 
GHz. We generated 100,000 impedance values. In addi- 
tion, to obtain the radiation impedance, we also used 
HFSS to simulate the case with radiation boundary con- 
ditions assigned to the sidewalls of the cavity. We find 
that the average value of the cavity reactance (which we 
predict to be the radiation reactance) has large system- 
atic fiuctuations. This is illustrated in Fig. O where we 
plot the median cavity reactance versus frequency. Here 
the median is taken with respect to the 100 locations of 
the perturbing disc. Also shown in Fig.Ois the radiation 
reactance Xfi{u!) = Im[Zii{u!)]. As can be seen the radi- 
ation reactance varies only slightly over the plotted fre- 
quency range, whereas the median cavity reactance has 




FIG. 2: (a) Median of ^ versus P/iV, compared with Eq. (UHl • 
(b) Histogram of approximation to P^{£,) (solid dots) in the 
TRS case compared with a Lorentzian distribution of unit 
width, (c) Same as (b) but for the TRSB case. 
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FIG. 3: Median cavity reactance averaged over 100 realization 
vs. frequencies ranged from 6.75GHz to 8.75GHz, compared 
with the corresponding radiation reactance Im[ZR{Lj)]. 



large frequency dependent fluctuations about this value. 
On the other hand, we note that over the range 6.75-8.75 
GHz, the average radiation reactance is 40.4 and the 
average of the median cavity reactances is 42.3r2. Thus 
over this frequency band, there is good agreement. The 
scale of the fluctuations in cavity reactance is on the order 
of 0.2 GHz, which is much larger than the average spacing 
between cavity resonances which is only 0.016GHz. Thus, 
these fluctuations are not associated with individual res- 
onances. Rather, the frequency scale of 0.2GHz suggests 
that they are multipath interference effects {L ^ lOOcm), 
which survive in the presence of the moveable conduct- 
ing disc. One possibility is that the fluctuations are the 
result of scars [s^l and this will be investigated in the 
future. The implication of Fig. 01 is that to obtain good 
agreement with the theory predicting a Lorentzian distri- 
bution, it may be necessary to average over a sufficiently 
large frequency interval. 

To test the Lorentzian prediction we normalize the 
cavity impedance using the radiation impedance as in 
Eq. 1251) and Eq. (|26|l . the normalized impedance values, 
^ = {Im[Z(A;)] — X R{k)]} / Rii{k) , are computed, and the 
resulting histogram approximations to P^iO is obtained. 
Fig. EJa) shows the result for the case where we have 
used data in the frequency range 6.75GHz to 8.75GHz 
(the range plotted in Fig. |3J). The histogram points are 
shown as dots, and the theorectical unit width Lorentzian 
is shown as a solid curve. Good agreement between the 
predicted Lorentzian and the data is seen. Figures 01 
(b)-(e) show similar plots obtained for smaller frequency 
range of width 0.5GHz: (b) 6.75 - 7.25 GHz, (c) 7.25 - 
7.75GHz, (d) 7.75 - 8.25 GHz, (e) 8.25 - 8.75 GHz. For 
these narrow freqency ranges, we see that Figs.01b) and 
^c) show good agreement with (|28|) . while, on the other 
hand, Figs.Qfd) and^fe) exhibit some differences. These 
are possibly associated with the variances in the median 
cavity reactance shown in Fig. 3 as the agreement with 
the Lorentzian prediction improves when averaging over 
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FIG. 4: Histogram approximation to P(^{^) from numerical 
data calculated using HFSS in different frequency ranges, (a) 
6.75 - 8.75 GHz, (b) 6.75 - 7.25 GHz, (c) 7.25 - 7.75GHz, (d) 
7.75 - 8.25 GHz, (e) 8.25 - 8.75 GHz. 
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FIG. 5: Schematic description of the two port model 
a large range of frequencies. 

C. Variation in Coupling 

In this section, we bolster our arguments connecting 
the radiation impedance and the normalization of the 
cavity impedance by showing that the relation is pre- 
served when the details of the coupling port are modi- 
fied. Let us consider a one-port coupling case in which 
the actual coupling is equivalent to the cascade of a loss- 
less two port and a "pre-impedance" Z seen at terminal 
2, as illustrated in Fig. |S1 

The impedance Z at terminal 2 then transforms to a 
new impedance Z' at terminal 1 of the two port according 
to 

Z =jAiiH , (30) 

JX22 + Z 

where jXij is now the purely imaginary 2 by 2 impedance 
matrix of the lossless two-port. We now ask how Z trans- 
forms to Z' when (a) Z is the complex impedance Zji 
corresponding to the radiation impedance into the cavity 
(i.e. the cavity boundaries are extended to infinity) and 
(b) Z = jX is an imaginary impedance corresponding to 
a lossless cavity, where X has a mean X and Lorentzian 
distributed fluctuation X. 

First considering case (a) the complex cavity 
impedance Zj^ = -\- jXji transforms to a complex 
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impedance Z't, = B!t, 



R'r = Rr 



jX'j^ where 

R\ + [X22 + XrY 



and 



X'ji — Xii — {X22 + Xr) 



X12X; 



21 



Rr 



{X2 



X, 



(31) 



(32) 



In case (b) we consider the transformation of the ran- 
dom variable X to a new random variable X' according 
to X' = Xii + Xi2i'2i/(-^22 + X). One can show that if 
X is Lorentzian distributed with mean Xr and width Rr 
then X' will be Lorentzian distributed with mean and 
the width i?^. Thus, the relation between the radiation 
impedance and the fluctuating cavity impedance is pre- 
served by the lossless two port. Accordingly, we reassert 
that this relation holds in general for coupling structures 
whose properties are not affected by the distant walls of 
the cavity. A treatment similar to that above has also 
been given by Brouwer |l5j in the context of scattering 
with a scattering matrix description of the connection 
between terminal 1 and 2. 

We now summarize the main ideas of this section. 
The normalized impedance of a lossless chaotic cavity 
with time-reversal symmetry has a universal distribution 
which is a Lorentzian. The width of the Lorentzian and 
the mean value of the impedance can be obtained by 
measuring the corresponding radiation impedance under 
the same coupling conditions. The physical interpreta- 
tion of this correspondence is as follows. In the radiation 
impedance, the imaginary part is determined by the near 
field, which is independent of cavity boundaries. On the 
other hand, the real part of the radiation impedance is 
related to the far field. In a closed, lossless cavity, the 
real part of the impedance vanishes. However, waves that 
are radiated into the cavity are reflected from the bound- 
aries eventually returning to the port and giving rise to 
fluctuation in the cavity reactance. 



IV. STATISTICS OF REFLECTION 
COEFFICIENT IN THE LOSSLESS CASE 

In the previous section, we obtained a universal 
Lorentzian distribution for the chaotic cavity impedance 
Z , after normalization by the radiation impedance. 



Z^j{XR + RRi), 



(33) 



where ^ is a zero mean, unit width Lorentzian random 
variable. We now consider the consequences for the re- 
flection coefficient. Suppose we can realize the perfect 
coupling condition, i.e. Rr ~ Zq, Xr = 0, in which 
the wave does not "feel" the transition from the cable to 
the cavity. In this case the cavity reflection coefficient 
becomes 



exp[-j(2tan"^^-h7r)]. 



(34) 



A standard Lorentzian distribution for ^ corresponds to 
a uniform distribution for tan~^ ^ from [— 7r/2, 7r/2], and 
thus to a refiection coefficient uniformly distributed on 
the unit circle. 

In the general case (i.e., non-perfect coupling), we in- 
troduce 7fl — Rr/Zq, 7x — Xr/Zq, and express S as 



S^e^^ = iZ+Zo)-'{Z-Zo) 



(35) 



We replace the Lorentzian random variable ^ by introduc- 
ing another random variable '0 via ^ = tan(?/'/2). Using 
this substitution, the Lorentzian distribution of ^ trans- 
lates to a distribution of -0 that is uniform in [0, 2n\. We 
then have from Eq. (|35|1 



J{<t>-'t>R) — 



1 



-jip' ' 



where the "free space reflection coefficient" pR 



PR. = \Pr\& 



IR + iix - 1 



IR + iix + 



1' 



(36) 



(37) 



is the complex reflection coefficient in the case in 
which the cavity impedance is set equal to the radia- 
tion impedance (^ = — j), and ip' = + "n + 4>r + 
2 tan"-^ ["fx /{iR + 1)] is a shifted version of i/i. Equations 
for the magnitude and phase of the free space reflection 
coefficient pR can be obtained from Eq. l|S7|l . Specifically, 



\PR\ 



ilR - 



l)^+7i 



(7i?. + l)'+7| 



2 ' 



and 



tan (l)R = 



27^ 



X 



-7i 



1 



(38) 



(39) 



Eq. (|36|l is essentially a statement of the Poisson kernel 
relation for a non-perfectly coupled one port cavity. 

To compute the probability distribution function for 
(/>, P^{(j)), we note that, since tp is uniformly distributed 
on any interval of 27r, we can just as well take tp' , which 
differs from -0 by a constant shift, to be uniformly dis- 
tributed. Consequently, we have 



27r ' d0 ' 
1 



(40) 



2nl + \pR\^-2\pR\cos{ 



Thus P(j){4>) is peaked at the angle (j)R corresponding to 
the phase angle of the free space reflection coefficient, 
with a degree of peaking that depends on \pr\, the mag- 
nitude of the free space reflection coefficient. 'Perfect 
matching' corresponds to 7^ = 1, jx = 0, and \pr\ — 0, 
in which case P,t,{<j)) is uniform. 
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We next consider the case of poor matching for which 
\Ppi\ = 1 and P^{(j)) is strongly peaked at (j)R- This be- 
havior can be understood in the context of the frequency 
dependence of the phase for a given reahzation. It fol- 
fows from (|35|l and (|27ll that the phase (j) decreases by 
27r as increases by the spacing between eigenvalues. If 
\pb\ = 1, then for most of the frequencies in this interval, 
the phase is near However, for the small range of 
frequencies near a resonance, the phase will jump by 27r 
as the resonance is passed. This indicates that the mode 
of the cavity is poorly coupled to the transmission line. 
In the case of good matching, = 0, all phases are 
equally likely indicating that, as a function of frequency, 
the rate of increase of phase is roughly constant. This 
implies that the resonances are broad, and the cavity is 
well coupled to the transmission line. 

In order to describe the different coupling strengths, we 
consider the parameter g originally introduced by Fyo- 
dorov and Sommers US : 



_ l + l(e^"^)P 

^ i-Ke^-^)r 

Evaluating (5*) using Eq. 



1 



\pr\ 



(41) 



(42) 



Thus, g has a minimum value of 1 in the perfectly 
matched case and is large if the matching is poor, ^ 
1. An analogous quantity is the voltage standing wave 
ratio on the transmission line when the cavity impedance 
is set equal to the radiation impedance. 



VSWR 



1 



\pr\ 



1 - \pr\ 



+ ^5 



1. 



(43) 



To test Eq. (|4(J|I . we compared its predictions for the 
phase distribution with direct numerical calculations ob- 
tained using HFSS (High Frequency Structure Simulator) 
for the case of the cavity and coupling detail as specified 
in Fig. ^ As compared to what was done for Fig. ^ we 
have narrowed the frequency range to 0.1 GHz bands for 
each realization in 1000 lO""* GHz steps centered at 7 
GHz, 7.5 GHz, 8 GHz, 8.5 GHz. Instead of calculating 
the radiation impedance for every frequency, we use the 
value of Zn at the middle frequency of the interval in 
calculating the values of 7^ and ^x- We present theo- 
retical phase density distribution functions together with 
numerical histogram results in Fig.|Bl The agreement be- 
tween the theory, Eq. 1)40(1 . and the numerical results is 
surprisingly good, especially considering the rather small 
(O.lGHz) frequency range used. 



V. EFFECT OF DISTRIBUTED LOSSES 

We now consider the effect of distributed losses in the 
cavity. By distributed losses, we mean losses that af- 
fect all modes in a frequency band equally (or at least 
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FIG. 6: Histogram of the reflection phase distribution for an 
HFSS calculation for the cavity in Fig. with center frequen- 
cies located at (a) 7GHz, (b) 7.5GHz, (c) 8GHz, (d) 8.5GHz, 
and with sweeping span equal to O.lGHz. Numerical data are 
compared with Eq. 14UII using parameters determined by Zr 
at the corresponding center frequencies. 



approximately so). For example, wall losses and losses 
from a lossy dielectric that fills the cavity are considered 
distributed. [For the case of losses due to conducting 
walls, the losses are approximately proportional to the 
surface resistivity, ~ ^//, and vary little in a frequency 
range A/ <C /. In addition, there will also be varia- 
tion of wall losses from mode to mode due to different 
eigenmode structural details. These modal fluctuations, 
however, are small when the modes are chaotic and the 
wavelength is short.] We use the random coupling model 
to construct a complex cavity impedance accounting for 
distributed losses in a manner analogous to the lossless 
case, Eq. (fTH|l . 



z^-^Va 

-TT ^ ' 



RR{kn)wl 



(44) 



where a represents the effect of losses. In particular, for 
loss due to wall absorption in a two-dimensional cavity, 
the value of a is equal to the ratio of the skin depth of 
the conductor to the height of the cavity; if the cavity 
contains a lossy dielectric, a is the loss tangent of the 
dielectric. The cavity quality factor is related to a by 
a — Q^^. This follows by noting that the real part of Z 
will have a Lorcntzian dependence on frequency (lu — kc) 
peaking at uj — A:„c with a full width at half maximum 
of Lua. 

The impedance Z will have a real part and an imag- 
inary part. We expect that, if k^a ^ A, corresponding 
to small losses, then the real part will be zero and the 
imaginary part will have an approximately Lorcntzian 
distribution. As losses are increased such that fc^cr ~ A 
(the imaginary part of the denominators in (|44() is of the 
order of eigenvalue spacing), the distributions of the real 
and imaginary part will change, reflecting that extremely 
large values of |Z| are no longer likely. In the high loss 
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limit, fc^cr ^ A, many terms in the sum contribute to the 
value of Z. In this case, we expect Z will approach the 
radiation impedance with small (Gaussian) fluctuations. 

In the Appendix we evaluate the mean and variance of 
the real and imaginary part of the complex impedance 
(|44|1 Z = R + jX. There it is shown that the mean is 
the radiation impedance Zj^ — Rn + jXf;, and the vari- 
ances of the real and imaginary parts are equal yar[i?] = 
V'ar[X]. In general, the distribution of R and X depends 
on the correlations between eigenvalues of k^. However, 
in the low damping limit, the correlations are unimpor- 
tant and we obtain 



Var[R] 



A 



27r fcV 



(45) 



for both the TRS and the TRSB cases. In the high damp- 
ing limit fe^cr 3> A, correlations are important and we 
obtain 



Var[R]^^^ 



for the TRS case 



for the TRSB case. 



(46) 



This is to be constrasted with the result one would obtain 
if correlations in the eigenvalue spacing were neglected; 
i.e., if the fc„ were assumed to be generated by adding 
independent spacings generated from the distributions 
(Uni and ((T7|l- In that case, using the method in the 
Appendix one obtains 



Var[R] 



A 

Rl A 



Z TT 

37r, 



for the TRS case 



(47) 



TT k^(7 16 



for the TRSB case. 



These results are larger than those in Eq. H46|l by 13.7% 
in the TRS case and 17.8% in the TRSB case, thus illus- 
trating the necessity of generating the fc^ using random 
matrix theory if accurate results are desired in the lossy 
case fc^cr > A. 

In a recent experimental paper the impedance 
statistics of a lossy TRS one-port microwave cavity were 
also considered. Their result is the same as (|44f) . One dif- 
ference is that they generate the realizations of fc^ solely 
by use of Eq. lfTC|l with the assumption that the eigen- 
value spacings are random independent variables. 

We now investigate a model, normalized impedance, 
applicable in the one-port case with loss, which is the 
generalization of Eq. (|27|) , 



C(^) 



N n 



kl 



(48) 



The normalized impedance C, will have a real part p > 
and an imaginary part ^, C = P + J?- We expect that if 
k^a <C 1, corresponding to small loss, then p = 0, and ^ 
will have an approximately Lorentzian distribution. 
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FIG. 7: (a) Histogram of the imaginary part of C,{(j) with 
different values of the damping obtained with method (ii); (b) 
Histogram of the real part of C,{a) with different dampings 
obtained with method (ii). (c) and (d) are histograms of 
the reactance and resistance from HFSS calculation with a 
lossy top and bottom plate, compared with histograms from 
Eq. H48|l computed as in (a) and (b) (dashed Une) and by 
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In analogy to Eq. H26|) we write for the cavity 
impedance 

Z = jXn + RrC, (49) 

and we use (|48|) to generate probability distribution func- 
tions for the real and imaginary part of C = p + j^. We 
first generate N values of Wn as independent Gaussian 
random variables of unit width (for this purpose we use 
a suitable random number generator) . We next generate 
N values of the normalized eigenvalues k'^. To do this we 
have utilized two methods: (i) an approximate method 
based on Eq. (jTHJ (for the TRS case) or Eq. in|) (for the 
TRSB case), and (ii) a method based on random matrix 
theory. We pick the value of relative to the spectrum 
fc^ such that the median of ^ is zero. 

For method (i) we independently generate N values 
of Sn using the distribution Ijltill or H17|l . We then ob- 
tain fc^ as fc^j = The main assumption of 
this method is that the spacings s„ can be usefully ap- 
proximated as uncorrelated. On the other hand, it is 
known from random matrix theory that the spacings are 
correlated over long distance (in n), and the thus the 
assumption of method (i) is questionable (compare (46) 
and (47)). This motivates our implementation of method 
(ii) (See also |3^). 

To implement method (ii) we generate an M x M ran- 
dom matrix with M large {M =1000) drawn from the ap- 
propriate ensemble (GOE or GUE) again using a random 
number generator. The width of the diagonal elements 
is taken to be unity. We then numerically determine the 
eigenvalues. The average spacing between eigenvalues of 
random matrices is not uniform. Rather, in the limit of 
large M, the eigenvalues A are distributed in the range 
— \/2M < A < \/2M, and the average spacing for eigen- 
values near an eigenvalue A is given by 

A(A) = 7r/V2A/- A2 (50) 

in both the TRS and TRSB cases. In order to gener- 
ate a sequence of eigenvalues with approximately uniform 
spacing we select out the middle 200 levels. We then nor- 
malize the eigenvalues by multiplying 1/A(0) to create a 
sequence of fc^ values with average spacing of unity. 

Histogram approximations to the GOE probability dis- 
tributions of Re[C] and /r7i[C] obtained by use of (|48|) and 
method (ii) are shown in Figs. 7(a) and 7(b). These 
were obtained using 30,000 random GOE matrix realiza- 
tions of 1000 by 1000 matrices and selecting the middle 
200 eigenvalues of each realization. The resulting graphs 
are shown for a range of damping values, fc^cr=0.01, 0.1, 
0.5, 1, 5 and 10. As seen in Fig. UJa), when k^a is in- 
creased, the distribution of f values becomes "squeezed" . 
Namely, the Lorentzian tail disappears and the fluctua- 
tions in ^ decrease. Eventually, when a enters the regime, 
1 <C k'^a <C N, the probability distribution function of 
^(cr) approaches a narrow Gaussian distribution centered 
at ^ = (recall that ^ = 0). As shown in Fig. [TJb), as 
a increases from zero, the distribution of the real part of 



C{<y) which, for cr = 0, is a delta function at zero, expands 
and shifts toward 1, becoming peaked around 1. In the 
very high damping case, both the real part and imagnary 
parts of ^, p and ^, will be Gaussian distributed with the 
mean value equal to 1 and respectively, and the same 
variance inversely proportional to the loss (as shown in 
the Appendix). As a consequence, the reflection coef- 
ficient \S\'^ in the high damping limit, is exponentially 
distributed. This result is consistent with the theoretical 
discussion given by [ssj . 

In general, the complex impedance distribution is not 
described using simple distributions such as Gaussian 
or Lorenzian. The distribution of the real part of the 
impedance has been studied in connection with the the- 
ory of mesoscopic systems and known as the "local den- 
sity of states" (LDOS). Through the supersymmetry ap- 
proach, Efetov obtained the probability density func- 
tion for the LDOS in systems without time reversal 
svmmetrv|2ll|. For chaotic systems with time reversal 
symmetry, the corresponding exact formula was derived 
in a form of multiple integral l34|. However the difli- 
culty to carry out the flve-fold integral makes it hard to 
interpret the formulus in Ref. |34| . Very recently, Fyo- 
dorov and Savin have proposed an interpolation formu- 
lus for the impedance distributions at arbitary values of 
damping parameter |35i| . The suggested formulas satisfy 
all the asymptotic behaviors in the physically interesting 
limiting cases, e.g. weak damping or very string damp- 
ing cases. Furthermore, these formulas seem to agree 
pretty well with the results of the numerical simulations, 
though the agreement in the intermediate damping case 
is not as good as in the limiting cases. In our paper, we 
still use the histograms generated from the Monte-Carlo 
simulations as a comparison to the HFSS data, however, 
we believe the formula presented by Fyodorov and Savin 
would be very helpful for the most of practical purposes 
of comparison. 

We noted that the variance of the real and imaginary 
parts of the complex impedance are equal. There is a 
more fundamental connection between these that is re- 
vealed by considering the reflection coefflcient in the per- 
fectly matched case, 

ae^"^ = (C-l)/(C + l), (51) 

where a and </> are random variables giving the magnitude 
and phase of the reflection coefficient. It can be argued 
|33| that (j) and a are independent and that (f> is uniformly 
distributed in [0, 27r]. The magnitude a is distributed on 
the interval [0, 1] with a density that depends on losses. 
A plot of the probability distribution for a taken from 
the data in Figs (7a) and (7b) is shown in Fig 8, for the 
damping values 0.1, 0.5, 1 and 5. 

We can express the actual complex reflection coefficient 
p in terms of the normalized reflection coefficient by first 
finding the normalized impedance from (|51|) . C = (1 -I- 
ae^'^)/{l — ae^'^) calculating the cavity impedance from 
(|49|l . and expressing the result in terms of the radiation 
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FIG. 8: Histogram of the magnitude of reflection coefficient 
in the Eq. I51II . a{<j), with different values of the damping. 



reflection coefBcient (|37ll . The result is 



PR 



1 + aeJ('^+^<^)p|j' 



(52) 



where tan(A(/)/2) = —Xji/{Rii + Zq) depends on system 
specific parameters. Since the angle <j) is uniformly dis- 
tributed, it can be shifted by Ai/) thus eliminating 
from the expression. Eq. 1)52(1 is then a restatement of 
the Poisson kernel in the single port case. 

The independence of a and (/) in Eq. H51() also 
guarantees the invariance of the distribution of cavity 
impedances when a lossless two port is added as in Sec 
III(C). In particular, the normalized cavity impedance C 
before the addition of the two port is given by 



c 



Z — jXfi 1 + ae- 



j't> 



Rr 



1 - ae-J"? 



(53) 



With the addition of the lossless two port as shown in 
the Fig. 5, impedances are transformed to Z' , X'^, and 



R'p^ such that 



c 



Z' - jX' 1 + ae^ 



1 - ae-J'C"*-"*-) ' 



(54) 



where 0c = (2/3 + tt) depends only on the properties of 
the two port and the cavity coupling port and the angle 
P satisfies 



cos/3 
sin/3 



Rr 



^Rj, + {Xn + Xh)^ 

(Xn + Xfl) 
v/i?| + {Xn + XrY ' 



(55) 



Since (j) is uniformly distributed, so is the difference — 
(j)c- Consequently, the normalized random variables C and 
C have identical statistical properties. 

A by-product of l(53|l is that we can easily prove that 
its real part p — {1 — a^)/(l + a"^ — 2a cose/)) and its 



imaginary part ^ = (2a sin 0)7(1 + — 2a cos 0) have 
the same variance and zero correlation. Since a and 
are independent, we can carry out the integration over 
the uniformly distributed (j) and obtain 



Var[p] = Var[^] = ( 



1 



1 — a- 



(56) 

where {■■)a denotes average over a. This property has 
been tested in microwave cavity experiments with excel- 
lent agreements |23|. For the high damping case, p — 1 
and ^ will become two independent Gaussian variables 
with zero mean and small but same variances. This 
yields an exponential distribution for the a^, which is 
consistent with the result obtained by Kogan based 
on the "maximum information entropy" principle. For 
the weakly absorpting case, Beenakker and Brouwer [s^ 
studied the distribution of a^ in the TRSB case through 
the time-delay matrix and obtained a generalized La- 
guerre ensemble. Howerver, for a TRS system with ar- 
bitary loss, there is no simple formula for the distribution 
of reflection coefficients. 

Using HESS, we simulate the lossy case by specifying 
the material on the top and bottom plates to be an imper- 
fect conductor with a bulk resitivity of 70 mVi-cm. In this 
case we can calculate a value of ct = 5/h = 0.002, where 
6 is the skin depth and h the cavity height. The corre- 
sponding parameter k^a is 0.5 at 7.75GHz. Histogram 
results for the normalized reactance (^) and resistance 
(p) fluctuations of Chfss = R'^{Zcav - jXr) = p + 
are plotted in Figs.^c) and^d) together with the his- 
tograms generated from Eq. (gHJ, and using spectra from 
the random matrices. As can be seen, the histograms 
from the HESS simulations match those of the model. 



VI. SUMMARY 

We have applied the concepts of wave chaos to the 
problem of characterizing the statistics of the impedance 
and scattering coefficient for irregular electromagnetic 
cavities with one port in the small wavelength regime. 
The coupling of energy in and out of the port in such 
cavities depends on both the geometry of the port and 
the geometry of the cavity. We found that these effects 
can approximately be separated. The geometry of the 
port is characterized by its radiation impedance which 
has both a real and imaginary part. This impedance de- 
scribes the port in the case in which the distant walls of 
the cavity are treated as perfect absorbers (or else are re- 
moved to infinity). The effects of geometry of the cavity 
can be treated in a statistical way using Random Marix 
Theory. The separation of the system specific aspects 
of the coupling and the universal aspects has previously 
been described using the Poisson kernel 0| . The relation 
of our approach to the Poisson kernel may be understood 
by comparing the equivalent relations (|49|l and (|52|l . To 
extract a universal quantity (^) from a set of impedance 



15 



values (Z) one must subtract the radiation reactance and 
normalize to the radiation resistance. To extract a uni- 
versal quantity (ae-''^) from a set of reflection values (p) 
we must solve the bilinear relation (|52|l for the magnitude 
a and phase (j) of the normalized reflection coefficient. 
The normalized impedance and scattering amplitude are 
related by ae^* = (C - 1)/(C + !)• If the radiation re- 
flection coefficient pu is known then (|52|l may be solved 
directly for normalized reflection coefficient, 



ae 



-.; A0 P ~ PR 



1 



PPr 



(57) 



The radiation reflection coefficient can be determined 
directly by measurement or by ensemble averaging. 
According to Ref. Q the average of p is equal to pn- 
This can be verified directly from H52|) by averaging over 
the uniformly distributed phase 0. Regardless of the 
value of a, one finds 



27r' 



PR- 



(58) 



Thus if enough appropriate, statistically independent re- 
alizations are available to computer the average of p, H58|l 
can be used to find the universal reflection amplitude. 

Consistent with previous results 0| our model pre- 
dicts that in the lossless case the impedance is Lorentzian 
distributed with a mean equal to the radiation reac- 
tance and a width equal to the radiation resistance. The 
Lorentzian prediction is tested by direct numerical solu- 
tion of Maxwell's equation for the cavity of Fig. The 
predictions are verified if an additional averaging over 
frequency is introduced. Effects of distributed loss and 
variation of coupling are also investigated and we have 
generated pdf's for the real and imaginary parts of the 
normalized impedance. In addition, we have calculated 
the mean and variance for these distributions and deter- 
mined the effect of correllations in the eigenfrequencies on 
the variances. The values of the variance depend on the 
degree of loss in the cavity and can be used to qunatify 
it. Finally, we have compared the predicted distributions 
of the normalized impedance with those obtained from a 
direct numerical simulation. 
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APPENDIX: VARIANCE OF CAVITY 
REACTANCE AND RESISTANCE IN THE 
LOSSY CASE. 

From Eq. H44() . we obtain the expression for the com- 
plex impedance in the single port case, 

. 1 ^ Mkl)RRikl)wl[kj+j{kl - P)] 

^ ' {k^-kiY + {klY ^ 



^Ria)+jX{a), 

(A.l) 

where A is the mean spacing (fc^ — X{a) and R{<t) 
are cavity reactance and resistance in the lossy case and 
= k^a. In this appendix, we are going to evaluate 
the mean and variance of X{a) and R{a) as well as their 
covariance. 

We first investigate the mean of R{a). We express the 
mean in terms of probability distribution function for the 
eigenvalues fc^. 



E[R{a)] 



1 



dkl . 



■ dk%Pj{kl 



^ RRA{wl,)kj 



(A.2) 



where 
{kl... 



Pj is the joint distribution of eigenlevels 
, k%) assuming they are unordered. Since the lev- 
els are not ordered, in each term of the sum, we can in- 
tegrate over all k"^ ^ /c^/, and obtain N identical terms. 
Thus, 



N f 

E[R{a)] = -J dkl,P^{kl,)RnMt 



k"^ 



(fc2 



U2 \2 



kj 

(A.3) 

where Pi{k^,) is distrubtion for a single level. Here we 
have introduced an integer N representing the total num- 
ber of levels. We will take the limit of — *■ oo. The 
single level probability distribution then satisfies by def- 
inition, 



Piikl) 



1 



NA{ki,y 



(A.4) 



We next assume that the radiation resistance RR{k^,) is 
relatively constant over the interval of fc^; values satis- 
fying — fc^, I < fc^ and we will move it outside the 
integral replacing it by Rii{k'^). Assuming that fc^ is not 
too large (fc^ <C fc^) we can take the end points at the 
integral to plus and minus infinity and evaluate Eq. I|A.3|) 
as 



E[R] 



Rr 

TT 



dx- 



1 



1 



i?i?(fc'), 



(A.5) 



where x — (fc^, — k'^)/k'^. Thus the expected value of the 
real part of cavity impedance is the radiation resistance 
independent of the amount of damping. This is some- 
what surprising since we have previously asserted that in 
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the lossless case, the cavity resistance is zero. The con- 
stancy of the expected resistance results from the reso- 
nant nature of the cavity impedance. When losses are 
small, fc^cr = fc^ ^ 1, for almost all frequencies the re- 
sistance is small. However, for the small set of the fre- 
quencies near a resonance the resistance is large. This 
is evident in the histograms of Fig. ^Tjp). The result is 
that small chance of a large resistance and large chance 
of small resistance combine to give an expected value re- 
sistance which is constant. 

In order to obtain the variance of i?(cr), we calculate 
the second moment of i?(cr), 

E[R{af] = f ... f dkl... dk%Pj{kl ...,k%) 



N 

E 



A^Rnikl,)Rn{kl,){wl,wl,)k^, 



= h+h- 

(A.6) 

Following the arguments advanced to calculate E[R{a)], 
we note that there will be N terms in the double sum for 
which kl^, — fc^, giving 



and N{N — 1) terms for which fc^, ^ fc^, giving 



N{N~1) / / dkl,dkl 



[{k'-ki,r + kj][{k^-ki,y + kj] 

(A.8) 

For the first integral we use (jA.4p for the single level dis- 
tribution function, and making the same approximation 
as before, we obtain 



T „2,,2, (^")A(fc^) 



(A.9) 



For the second integral we need to introduce the two level 
distribution function. For the spectra that we consider, 
the two level distribution has the form 



p,(ki,,ki,) = i^ni 



(A.IO) 



Here the function g describes the correllations between 
two energy levels. For uncorrellated levels giving a Pois- 
son distribution of spacings we have 5 = 0. In the 
presence of level repulsion we expect 5(0) = 1 with 
{1 — g) (X |/c^, — k'^,\^ for small spacing, and /3 = 1 for 
TRS and /3 = 2 for TRSB systems. As \kl, -kl,,\^ 00, 
g ^ indicating loss of correllation for two widely sep- 
arated levels. The function g will be different for spec- 
tra produced by random matrices and spectra genergated 
from sequences of independent spacings. Expressions of 
g for the spectra of random matrices can be found in 



the book by Mehta ([13, Ch. 5 & 6). We will derive 
the expression for g for spectra generated by sequences 
of independent spacings later in this appendix. 

Based on expression (jA.10|l and the usual assumptions 
on the slow variations of Rn and A with eigenvalue fc^, 
we obtain 



h = {E[R]f-Ig 



(A.ll) 



where the first term comes from the 1 in lA.lOl and the 
second term comes from the correllation function g 



RRik^){w 



2^/,„2\2 /.oo 



ki A+d^kif 



The variance of R is thus given by 
Var[R] = E[R]^ - E[R^] 



(A.12) 



dfc2 2g(|fc2|) 



(A.13) 

Note, since w is a Gaussian random variable with zero 
mean and unit variance, (w^) = 1 and (w^) = 3. 

Equation (jA.lSp shows that the variances of R depends 
on fc^/A, the ratio of the damping width to the mean 
spacing of eigenvalues. In the low damping case, A:^/A <C 
1, the integrand in (|A.13|I is dominated by the values of 
|fc2| < A and we replace g by its value .g(0). Doing the 
integral we find 



Var[R] = R\ 



A 



fr2 



27r 



-gmw^ 



(A. 14) 



Since the damping is small, the first term dominates and 
the variance is independent of the eigenvalue correllation 
function. This is consistent with our previous findings 
that the eigenvalue statistics did not affect the distribu- 
tion of reactance values. 

In the high damping limit, A:^ > A, the integral in 
(jA.13|l is dominated by k"^ values of order A and we have. 



I 



(A.15) 



The variance decreases as damping increases with a coef- 
ficient that depends on the correllation function. Physi- 
cally the correllations are important because in the high 
damping case a relatively large number of terms in the 
sum (lA.lll contribute to the impedance and the sum is 
sensitive to correllations in these terms. 

The integral of the correllation function can be eval- 
uated for different spectra. For spectra generated from 
random matrices, we have f|l3j|. Ch.6) 

= f{sf - r^nis)] (A.16) 



for TRS matrices and 



9is) = f{sf 



(A.17) 
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for TRSB matrices, where /(s) = sin(7rs)/(7rs). In both 
cases, we find 



dsg{s) 



1 



(A.18) 



However, to consider the TRSB case we need to repeat 
the calculation including complex values of the Gaussian 
variable w. The result is 



°° dlf_ 
~^ 



(A.19) 



For spectra generated by sequences of independent spac- 
ing distributions we will show 



^^^9{\k'\) = l-\{s'), (A.20) 



where (s^) is the expected value for the normalized near- 
est neighbor spacing squared. Using and ((T7|l . this 
gives 



^9{\k'\) 



1 
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for TRS, 
for TRSB. 



(A.21) 



Note also that ljA.20p gives the required value of zero for 
Poisson spacing distributions, where (s^) = 2. 

We can evaluate the expected value of the reactance 
and its variance, as well as the covariance of reactance 
and resistance, using the same approach. We find the 
expected value of reactance is given by the radiation re- 
actance. 



E[X]=XR{k^). 



(A.22) 



The variance of the reactance is equal to that of the re- 
sistance (TOal the covariance between them is zero. 

We now derive the g-integral (|A.20|) for spectra gener- 
ated from independent spacings. We introduce a condi- 
tional distribution Pm{s) that is the probability density 
that the m}^ eigenvalue is in the range [s, s + ds] given 
that eigenvalue m = 0, is at zero. For convenience, here s 
is the normalized spacing with unit mean. When m = 1, 



Pi(s) is the spacing distribution p{s). Thus, 1 — g{s) 
stands for the probability that there exists an eigenlevel 
at [s, s-fds] given one level located at 0. This equality 
can be expressed as the summation of P.„i{s), 



(A.23) 



Pm{s) can be evaluated assuming the spacings are inde- 
pendent, 



CO in IIU 

1 - = E [ / n ds^Pl{s.)5{s - s,)]. (A.24) 

ni—1 i—1 i— 1 

We Laplace transformc both sides of Eq. 1A.24I) , and ob- 
tain 

(A.25) 

To evaluate dsg{s), we take the limit of t ^ 0. The 
transform -Pi(t) can be expressed in terms of the mo- 
ments of Pi{s), 



Pi(t) = / e-'^Pi{s)ds, 



foo 2 2 

/ (l-,sr+i— (A.26) 
Jo ^ 



= l-r{s) + -r{s') 



Thus, we can evaluate the integration of g{s) to be 

I dsg{s) — lim / dse^"^^ g{s) 
Jo Jo 



= lim[ 



A(r) 



oV l-Pi(r) 



(A.27) 



which is Eq. (|A.20|I 
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